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Abstract. The strong gravitational field of neutron stars in the brany universe 
could be described by spherically symmetric solutions with a metric in the exterior 
to the brany stars being of the Reissner-Nordstrom type containing a brany tidal 
charge representing the tidal effect of the bulk spacetime onto the star structure. We 
investigate the role of the tidal charge in orbital models of high-frequency quasiperiodic 
oscillations (QPOs) observed in neutron star binary systems. We focus on the 
relativistic precession model. We give the radial profiles of frequencies of the Keplerian 
(vertical) and radial epicyclic oscillations. We show how the standard relativistic 
precession model modified by the tidal charge fits the observational data, giving 
estimates of the allowed values of the tidal charge and the brane tension based on the 
processes going in the vicinity of neutron stars. We compare the strong field regime 
restrictions with those given in the weak-field limit of solar system experiments. 
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1. Introduction 

The string and M-theories describing gravity as higher-dimensional interaction 
appearing to be effectively 4D at low energies inspired braneworld models with 
the observable universe being a 3-brane (domain wall) to which the standard non- 
gravitational matter fields are confined; on the other hand, gravity enters the extra 
spatial dimensions allowed to be much larger than Zp ~ 10~ 33 cm pQ. Gravity can be 
localized near the brane even with an infinite size extra dimension with the warped 
spacetime satisfying the 5D Einstein equations with a negative cosmo logical constant, 
as shown by Randall and Sundrum [2]. Consequently, an arbitrary energy- momentum 
tensor could be allowed on the brane with effective 4D Einstein equations [3]. 

The Randall-Sundrum model gives standard 4D Einstein equations in the low- 
energy limit, but significant deviations occur at very high energies, e.g., in the vicinity 
of compact objects, such as black holes jl] or neutron stars [5]. There are high-energy 
effects of local origin influencing pressure of matter and non-local effects of "back- 
reaction" origin arising from the influence of the Weyl curvature of the bulk space on 
the brane. The combination of the high-energy (local) and bulk stress (non-local) effects 
alters significantly the matching problems on the brane in comparison with the standard 
4D gravity |6j, and the stresses induced by the bulk gravity imply that the matching 
conditions do not have unique solution on the brane [6]. The 5D Weyl tensor is needed 
as a minimum condition for uniqueness, but no exact 5D solution is known recently [5]. 

Assuming a spherically symmetric metric induced on the brane, solutions of the 
effective 4D Einstein equations could be found relatively easily. The vacuum (in both 
the bulk and brane) solutions represent black holes (naked singularities) of the Reissner- 
Nordstrom (RN) type endowed with a brany tidal charge parameter b, describing the 
bulk tidal stress effect on the black hole structure [I] . In the brany RN metric, the tidal 
charge b stands instead of the charge parameter Q 2 of the standard RN spacetimes [7]; b 
can be both positive and negative, but there are indications that b < should properly 
represent the effects of the bulk space Weyl tensor [H 0, [6j 

In the case of spherically symmetric neutron stars (or quark stars) simplified by the 
uniform energy density profile, the local high-energy and non-local bulk graviton stress 
effects make the matching conditions unambitious, since we do not know the 5D Weyl 
tensor. Note that the uniform energy density distribution model of neutron stars could, 
in principle, be used in order to roughly determine the neutron star spacetime structure, 
since in realistic models with a variety of realistic equations of state, the energy density 
profile happens to be nearly uniform in most of the neutron star interior jTQl ITT] . 

Two different exact exterior solutions related to the uniform density interior are 
known, both having asymptotically Schwarzschild character [5]. The first one is of the 
known RN type with the tidal charge and mass parameters being influenced by the 

| The stationary and axially symmetric vacuum solutions describing rotating black holes (naked 
singularities) are found to be of the Kerr-Newman type endowed again with the tidal charge reflecting 
the effect of the bulk stresses [9]. 



QPOs in brany neutron stars 



3 



energy density of the star and the brane tension. The second one is more complicated, 
with metric components being explicitly dependent on the energy density and brane 
tension [5] . Quite recently, the solution of the RN type was used in the weak-field limit 
in order to obtain restrictions on the brany parameter and the brane tension from the 
standard tests on the Mercury perihelion precession, deflection of light near the edge of 
Sun and the radar echo experiments [12J. Here, we also restrict our investigations to 
the metric of the RN type in order to obtain an idea of the tidal charge relevance in a 
strong gravitational field near neutron stars (or black holes) and to find restrictions on 
the allowed values of the brany parameter. The choice of the RN type of the uniform 
star brany solutions is supported by a recent study of static exteriors of non-static brany 
stars showing that spherical non-rotating star solutions admit only Schwarzschild or RN 
type exterior [13]. Since the junction conditions of the internal and external neutron 
star spacetimes relate the tidal charge to the brane tension and the energy density (mass 
and radius) of the neutron star, we can, in principle, also put restrictions on the allowed 
values of the brane tension [5]. 

There is a variety of papers related to both weak-field [H] and strong-field 
[T5| [IHl [T7j optical lensing phenomena in the braneworld black hole spacetimes. The 
optical phenomena connected to the accretion disc radiation could also reflect in a 
relevant way the influence of the tidal charge [18j [19] . Here we focus our attention on 
the strong-field phenomena closely related to the perihelion precession effect discussed 
in the weak-field limit [12] . 

High-frequency quasiperiodic oscillations (QPOs) observed in Galactic low-mass 
binary systems, namely, in both the black hole [20l |2T| [22l [23] and neutron star binary 
systems [241 EHJ EHJ [27] , or in extragalactic black hole sources [28], [29], EDI EI] , represent 
such an interesting phenomenon in the strong field regime enabling estimates of the 
tidal charge on the basis of orbital models. It is widely accepted that the high-frequency 
QPOs are related to the orbital motion in the inner part of an accretion disc around a 
compact object. There are strong indications that some resonant phenomena could be 
relevant in both black hole systems where twin peak QPOs with fixed frequencies and 
their ratios are observed [23], and neutron star systems with wide spread of frequencies 
and their ratios where the twin peak QPOs are cumulated around small number ratios 
[321 1.33J. The hypothetical resonances of the orbital epicyclic oscillations representing 
plausible models of QPOs were recently studied for brany Kerr black holes [31]. Here, 
we study these models as related to the brany neutron stars described approximately 
by the RN spacetimes, i.e., neglecting the effects related to rotation of neutron stars. 
We focus our attention to the relativistic precession model [35] that fits qualitatively 
well in the wide spread data given by observations of a variety of neutron star binary 
systems [2U [251 I2S1 [27]. However, recently it has been discussed in the case of the 
atoll source 4U 1636—53 [32J E3J [36] and some other atoll sources [37] that the resonant 
phenomena should be relevant even in the framework of the relativistic precession (RP) 
model. We, thus, present a first study of strong-field phenomena for brany neutron stars 
that could be compared with the weak-field results (perihelion precession, radar echo 
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delay and light deflection) based on the RN solutions used to describe approximately the 
spacetime in the vicinity of Sun [12]. In order to obtain rough restrictions on both the 
brany tidal charge and brane tension, we use an ensemble of typical neutron star atoll 
sources and the special cases of the Z-sources Sco X-l and Circinus X-l that demonstrate 
a wide variety of QPOs. 

In section [2J the effective 4D Einstein equations are presented, and the braneworld 
neutron star spacetime of the RN type (approximated by the uniform energy density 
distribution |5J) is briefly described. In section [3], the circular geodesies are discussed. 
In section HI frequencies of the orbital epicyclic oscillations are given; their radial profiles 
are discussed, and the resonance conditions are briefly described for some frequency 
ratio functions modeling twin peak QPOs. In section [5l the RP model is used to fit 
data observed in the representative ensemble of neutron star X-ray binaries exhibiting 
QPOs, and the restrictions on the tidal charge and brane tension are given. In section El 
concluding remarks are presented. 



2. Gravitational field equations on the brane and the exterior spacetime to 
the uniform density configurations 

In the 5D warped space models involving a non-compact extra dimension, the 
gravitational field equations in the bulk can be expressed in the form [31 H] 

G AB = k 2 [-Ag AB + 5( X )(-Xg AB + T AB )], (1) 

where the fundamental 5D Planck mass Mp enters via k 2 = 87r/M|, A is the brane 
tension, and A is the bulk cosmological constant. 

The effective gravitational field equations induced on the brane are determined by 
the bulk field equation (TjQ), the Gauss-Codazzi equations and the generalized matching 
Israel conditions with Z 2 symmetry. They can be expressed as modified standard 
Einstein's equations containing additional terms reflecting bulk effects onto the brane 
®: 

= -Kg^ + k% u + k 2 S^ - E^, (2) 

where k 2 = 8ir/Mp, with Mp being the brany Planck mass. The relations of the energy 
scales and cosmological constants are given in the form 



Mp = i/-H^)Mp; A = ^ 
^\VXJ ' M| 




3M| 



(3) 



Requirement of zero cosmological constant on the brane (A = 0) implies the relation 
A=-«£. (4) 

Local bulk effects on the matter are determined by the "squared energy-momentum" 
tensor arising from the extrinsic curvature term in the projected Einstein tensor: 

5^ = yTT» v - l -T«T va + 1^ (3T^T aP - T 2 ) , (5) 
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where T^ v is the energy-momentum tensor on the brane. The non-local bulk effects are 
given by the tensor E^ v representing the bulk Weyl tensor, Cacbd, projected onto the 
brane, whereas 

E AB = C A cbd n c n D . (6) 

There are two known solutions of the effective 4D Einstein equations corresponding 
to the brany uniform density configurations that could roughly represent the spacetime 
exterior to neutron (quark) stars [5]. Here, we focus our attention on the solution of the 
Reissner-Nordstrom type that is formally identical to the spherically symmetric black 
hole solution; however, its parameters are related to the brane tension A and energy 
density p of the internal configuration due to the matching conditions of the internal 
and external spacetimes. The external spacetime is given in the standard Schwarzschild 
coordinates by the line element: 

ds 2 = - A(r)c 2 dt 2 + A- 1 (r)dr 2 + r 2 d^ 2 + r 2 sin 2 ^d0 2 , (7) 

Ar=l-- + i; ^g = — , (8) 

where M is the gravitational mass, and B is the brany tidal charge representing the bulk 
non-local effects on the 4D spacetime structure. The matching conditions imply [5] 

*=-^*ra (9) 



c 2 VA 

m = m(i-£), (io) 

where R is the internal configuration radius, and the effective mass is given by 

R 

M = An J p cS (r)r 2 dr, (11) 
o 

with the effective density given by 

+ + < 12 > 

The non-local bulk gravitational effects arising from the bulk Weyl tensor are represented 
by the "dark energy" density U. The standard general relativistic equations are regained 
in the limit of A -1 — > 0. 

The internal solution puts two important restrictions on the brane tension A [5]. 
The first one is general for all uniform stars and reads 

Aa ( fi G 2GM/cO P - <13) 
implying R > 2GM/c 2 , i.e., the Schwarzschild radius is still the relevant limit as in 
general relativity. 

The second one represents an upper limit on compactness of the star following from 
the requirement that pressure must be finite inside the star. It takes the form [5] 



GM/c 2 4 



1 + 5p/4A 
R ~ 9 L(l + p/A) 2 J 



< 



(14) 



QPOs in brany neutron stars 



6 



that implies the general relativistic limit (GM/c 2 ) /R < 4/9 for A -1 — > 0. We can see 
that the braneworld high-energy corrections reduce the compactness limit of the star 
and the lowest order correction reads 



3. Circular geodesies of the external spacetime 

Considering vacuum spacetimes, the event horizons of the brany RN metric are 
determined by the condition A(r) = 0. Using geometric units (c = G = 1), the radius 
of the outer event horizon is given by the relation 



The horizon structure depends on the sign of the tidal charge B. We see that, in contrast 
to its positive values, the negative tidal charge tends to increase the horizon radius (see, 
e.g., figured]). 

It is clear that the positive tidal charge acts to weaken the gravitational field, and we 
have the same horizon structure as in the usual Reissner-Nordstrom solution. However, 
new interesting features arise for the negative tidal charge. For B < 0, it follows from 
equation (fl6l) that the horizon radius [4J 



such a situation is not allowed in the framework of general relativity. 

The event horizon does exist provided that M 2 > B, where the equality corresponds 
to the extreme case. For B > M 2 we obtain a naked singularity spacetime. Clearly, 
such a situation is impossible for negative values of the brany tidal charge B. Here, we 
are confronted with another situation, since the vacuum spacetime is limited to r > R 
and at r = R it has to be matched to the internal uniform density spacetime. Clearly, 
the radius R must satisfy the compactness condition (fl4|) . 

Note that values of the brany tidal charge B are not restricted for both positive 
and negative values; therefore, the external spacetimes with B > M 2 (corresponding to 
the "naked singularity" case) have to be considered too. 

In the analysis of the epicyclic frequency profiles, it is useful to relate the profiles to 
the photon circular geodesic and innermost stable circular geodesic radii (or innermost 
bound circular geodesic radii in the case of thick discs not considered here), which 
are relevant in the discussion of properties of the accretion discs and their oscillations. 
Therefore, we put the limiting radii in an appropriate form. 

Introducing the dimensionless radial coordinate 




(15) 



M + VM 2 - B . 



(16) 



r+ > 2M; 




r 

x = — 



(18) 



and a dimensionless brany parameter 



B 
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the metric coefficient takes the form 

A(x,b) = l-- + \. (20) 

x x 2 

The motion of a test particle with mass m is given by the standard geodesic equation 
DC/^/dr = and the normalization condition U^U^ 1 = —m 2 . The motion is restricted 
to the central planes; for one particle the plane can be chosen to be the equatorial plane 
(6 = 7r/2 = const). Then, the relevant radial equation of motion takes the form 

(hV = E 2 -V? s (x,L,b), (21) 



dr 

where the effective potential 

V 2 s (x, L, b) = A(x, b) (m 2 + ^ . (22) 

There are two constants of motion related to the spacetime symmetries: E/m is 
the specific energy; L/m is the specific angular momentum. The circular geodesic 
orbits (determined by dV e fi/dx = 0) have the specific energy and the specific angular 
momentum given by 

E = x *-2x + b 

m x {x 2 -3x + 2b) 1/2 ' 



L = ± (24) 

m (x 2 -3a; + 26) 1/2 



From equations ([23]) and ff24|) we can see that the circular orbits can exist from 
infinity down to the radius of the limiting circular photon orbit, determined by the 
condition 

x 2 - 3x + 2b = 0, (25) 

where E/m — > oo and L/m — > ±oo, but the impact parameter I = L/E remains finite. 
For "naked singularity" RN spacetimes, the restriction 

x > b (26) 

can be relevant. 

The radius of the marginally bound orbits with E 2 = m 2 is given by the largest 
root of the equation 

x (Ax - x 2 - 46) + b 2 = 0, (27) 
while the stable circular orbits are restricted by the condition 

x (Qx - x 2 - 9b) + Ab 2 < 0. (28) 
The outer event horizon x^(b) is implicitly determined by the relation 

b = b h (x) =x{2-x), (29) 
the radius of the photon circular orbit x p h(6) is given by 

b = b ph (x) = ^(3-x), (30) 
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the radius of the marginally bound orbit x m b(6) is given by 

b = 6 mb (x) = x (2 T y/x) , (31) 
and the radius of the marginally stable orbit x ms (6) is given by 



b = 6 ms (x) ~ | (9 T Vl6x-15) . (32) 

o 

The upper sign in (I3"T|) and (I32p is relevant for the black hole spacetimes in the 
region above the outer horizon, while both signs are relevant for the naked singularity 
spacetimes. 

Now, we have to put the limits on the applicability of the vacuum spacetime using 
the limit condition on the compactness of the uniform internal configuration given by 
the relation (|14p . We shall give the limits in terms of the dimensionless brany tidal 
charge b that can be expressed in the form 

&=-3-(T)=-3x(f). (33) 



r G VA/ VA< 

Clearly, negative (positive) values of b correspond to positive (negative) values of the 
brany tension A. 

The compactness relation (TI4~|) then reads 



and in the linear form 



"77 ) — 7T~~r ■ (35) 



RJ ~ 9-6 

The equality in relation (l34"j) determines the critical radius Rq representing the limit on 
the radius of neutron stars. We have to restrict our studies to the region r > R c . The 
radii Xq = Rc/rc are illustrated in figured! and are given (in an implicit form) by the 
relation 

b = b c (x) = - (I8x - 5x 2 =f x 3/ V25x- 36) . (36) 
6 

The functions bh(x), b p h(x), b mh (x), b ms (x) and be (x) are illustrated in figured! It is 
evident that the positive tidal charge will play the same role in its effect on the circular 
orbits as the electric charge in the Reissner-Nordstrom spacetime - the radius of the 
circular photon orbit, as well as the radii of the innermost bound and the innermost 
stable circular orbits shift inwards as the positive tidal charge increases. For the negative 
tidal charge the limiting photon orbit, the innermost bound and the innermost stable 
circular orbits shift outwards as the absolute value of b increases jU [9]. Of special 
interest is the behaviour of bc{x) putting limits on validity of the internal uniform 
density configuration. For 6 = 0, the limit is Xq = 9/4, and it falls with 6 > growing 
up to 6 = 2.6, where it starts to grow again. On the other hand, for 6 < the critical 
radius X c slowly grows with 6 descending to higher negative values, but for 6 = —1.5, 
there is X c = x^, and the function bc(x) loses its relevance leaving X > x^ as the only 
basic restriction on the surface radius of the neutron star. 
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The function b^x) has a local maximum at x = 1 where 6=1, corresponding 
to the extreme case. This local maximum coincides with the local minimum of the 
function 6 ms (x) at x = 1. The function 6 p h(x) has a local maximum at x = 3/2 for 
6 = 9/8 corresponding to the naked singularity spacetimes. The local maximum of 
the function 6 m b(x) is located at x = 16/9 for 6 = 32/27; the function 6 ms (a;) has a 
local maximum situated at x = 5/2 for 6 = 5/4 (this local maxima appear again only 
for naked singularities); the function b c (x) crosses bh(x) at 6 = —1.5 (see figure [p. 
Clearly, due to the compactness limit bc(x), the photon circular orbits are relevant in 
the spacetimes with 6 < 1.017 only. Further, we immediately see that for 6 > 0, the 
marginally stable orbits are well defined in the black hole type spacetimes (6 < 1) and 
in naked singularity type spacetimes with 6 < 5/4 and 6 > 9/4. In the spacetimes with 
5/4 < 6 < 9/4, the stable orbits can be located at all x > Xc, since x ms < Xc (and 
x ms = Xq for 6 = 9/4) in such vacuum RN spacetimes. 



4. Epicyclic oscillations of Keplerian motion 

It is well known that for oscillations of both thin Keplerian [381 [39] and toroidal discs 
[40] around neutron stars (black holes) the orbital Keplerian frequency z/r and the 
related radial and vertical epicyclic frequencies v T and ug of geodetical quasicircular 
motion are relevant and observable directly or through some combinational frequencies 
[231 HU H2], H3] . Of course, for extended tori, the eigenfrequencies of their oscillations are 
shifted from the epicyclic frequencies in dependence on the thickness of the torus [44~| 145] . 
Similarly, due to non-linear resonant phenomena, the oscillatory eigenfrequencies could 
be shifted from the values corresponding to the resonant geodetical epicyclic frequencies 
in dependence on the oscillatory amplitude [461 H7]. It is expected that shift of this kind 
is observed in neutron star systems [481 149] . while in microquasars, i.e., binary black 
hole systems, the observed frequency scatter is negligible and the geodetical epicyclic 
frequencies should be relevant. 

In the case of the external braneworld neutron star spacetime of the RN type with 
the brany tidal charge 6, the formulae of the test particle geodetical circular motion and 
its epicyclic oscillations, obtained in [50], could be directly applied (putting a = and 
Q 2 = 6). We can write the following relations for the orbital and epicyclic frequencies: 

1 /GM\ 1/2 1 / 6\ 1/2 Wi_^Y /2 

UK ~ 2tt V r 3 G J x 3 /2 { x) 2tt GM x^ \ x J ' ( j 

v l = a^K; ( 38 ) 
= a x v\, (39) 

where the dimensionless quantities determining the epicyclic frequencies are given by 

a e = 1, (40) 

(, &VV-, 6 9 6 46 2 \ 
«r= l-- 1 - - + -5 - -ST , (41) 



XI \ XX 1 x 3 



so that vk(x, 6) = vg(x, 6) due to the spherical symmetry of the spacetime. 




Figure 1. (a) The functions bh (black solid line), b p h (dashed-dotted line), 6 m b (dashed 
line) and b ms (thin black solid line) that implicitly determine the radius of the outer 
event horizon, the limiting photon orbit and the marginally bound and marginally 
stable circular orbits of a braneworld RN type spacetime. The grey solid line represents 
the function be implicitly determining the critical radius Xq of neutron stars. The 
thick black solid lines represent the functions Bk{x = Xk) and (x = X T ), implicitly 
determining the locations Xk and X r of the Kcplerian and radial epicyclic frequency 
local extrema. Thick grey dotted line represents the function b = x determining the 
other limit on the reality of circular orbits. The function 6 sr that implicitly determines 
the strong resonant radius x sr , where v^(x,b) = u T (x,b) (a T (x,b) = 1) is denoted 
by thick grey dashed line. The function & 3:2 (black dotted line) represents the radii, 
where the relativistic precession resonance i>k '■ (^k — v T ) =3:2 occurs, (b) The 
detailed figure introduces the relevant values of the brany tidal charge, separating the 
external spacetimes of different properties. 
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In the field of brany RN black holes (neutron stars) with b < 1, there is (see figure [2]) 



However, this statement is not generally correct in the case of brany RN naked 
singularities (neutron stars) with b > 1. For b > 1.42, the case u K (x,b) < u T (x,b) 
is even possible in some region of the external (vacuum) spacetime. 

Now we can discuss the behaviour of the fundamental orbital frequencies for 
Keplerian motion in the field of both brany RN black hole and brany RN naked 
singularity spacetimes, i.e., we consider external neutron star spacetimes with tidal 
charge b < 1 and b > 1. (Note that b < always corresponds to the black hole 
spacetime behaviour.) 

We express the frequency as v [Hz] 10 M©/M in every quantitative plot of frequency 
dependence on radial coordinate x, i.e., displayed value is the frequency relevant for a 
central object with a mass of 10 M©, which could be simply rescaled for another mass 
just by dividing the displayed value by the respective mass in units of ten solar masses. 

4-1. Properties of the Keplerian and radial epicyclic frequency 

First, it is important to find the range of relevance for the functions z/k(x,6) and 
u T (x,b) above the critical radius Rq of neutron stars. For completeness we consider 
their behaviour above the event horizon Xh for black holes (neutron stars), and above 
the ring singularity located at x = for naked singularities; recall that in the external 
spacetimes with b < —1.5, the restriction by x > Xc(b) is irrelevant, but the horizon 
condition x > %(&) keeps its relevance. 

The circular geodesies can exist for x > x p h(b), where x p h(b) is implicitly determined 
by the equation ( 1251) . Stable circular geodesies, relevant for the Keplerian, thin accretion 
discs, exist for x > x ms (b), where x ms (b) is determined (in an implicit form) by the 
relation (12"8]1 . which coincides with the condition 



For toroidal, thick accretion discs the unstable circular geodesies can be relevant in the 
range x m b < x in < x < x ms , being stabilized by pressure gradients in the tori [HU 152] . 
The radius of the marginally bound circular geodesic x m b, implicitly determined by the 
equation (!27l) . is the lower limit for the inner edge of thick discs (53j [54]. 

Clearly, the Keplerian orbital frequency is well defined up to x = x p ^. However, v r 
is well defined, if a T > 0, i.e., at x > x ms , and u T = at a; ms . From figure [21 we can 
conclude that not only the radial epicyclic frequency but even the Keplerian frequency 
can have a maximum located above the outer event black hole horizon. In the following 
subsection, we will discuss whether the maximum of vk{x, b) could be located above the 
marginally stable or the limiting photon circular orbit. 



z/ K (x, b) > v T (x, b). 



(42) 



a T (x, b) = 0. 



(43) 
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x 



Figure 2. The behaviour of the Keplerian frequency uk and radial epicyclic frequency 
v t in the field of braneworld RN black hole (neutron star) spacetimes with various 
values of the tidal charge parameter b. The curves are spaced by 0.2 in b, and they 
are plotted from the outer event black hole horizon x^. The dashed lines represent 
Schwarzschild spacetime with zero tidal charge. 



4-1.1. Local extrema of the Keplerian frequency. We denote by Ak the local extrema 
of the Keplerian frequency vk(x, b) given by the condition 

From ( 1371) . we find that after introducing ' as d/dx the corresponding derivative is given 
by 

1 IGM Ab~3x Ab-3x 

and relation ( 144"1) implies that the Keplerian frequency has a local extremum located at 

X K = (46) 
The second derivative at x = Ak, 

„_ 1 [GMJU/3_ 

% ""^V^i^v!' (47) 

is always negative; thus, the Keplerian frequency has a local maximum at x = Ak- 

For 6 = 15/16, the maximum of the Keplerian frequency radial profile is situated 
exactly at the radius coinciding with the radius of the outer event horizon Ak = Xh = 5/4 
(see figures 1(b) and [3]). For b = 9/8, the maximum of the Keplerian frequency is 
located at the circular photon orbit = x p h = 3/2 (this case corresponds to the naked 
singularity spacetime). 
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We can conclude that for brany parameter from the interval 

15/16 < b < 9/8 (48) 

the Keplerian frequency radial profile has its maximum located between the outer event 
horizon and the circular photon orbit: 

x h < X K < x ph , (49) 

i.e., it is physically irrelevant there. 

Generally, the maximum of the Keplerian frequency is physically irrelevant for all 
the brany RN neutron stars (of the black hole type with b < 1), since it could never be 
located above the circular photon orbit x p h (and the marginally stable orbit x ms ). For 
spacetimes with b < 1, the Keplerian frequency is a monotonically decreasing function 
of radius in astrophysically relevant radii above the photon orbit, x > x p h- The situation 



is clearly illustrated in figure 1(b). For the RN naked singularity type spacetimes with 
b < 1.381 the extrema of v& must be hidden under the neutron star surface. The local 
extrema of z/k can appear in the external spacetime to the internal uniform density 
configuration for b > 1.381 (see also figure [3]). 

4-1.2. Local extrema of the radial epicyclic frequency. The local extrema of the radial 
epicyclic frequency X r are given by the condition 

Using ([MP , the derivative reads 

a' 



v[ = ^ K + ^-z/ K J , (51) 
, 86 5 1 — x 1 „ . 

-^ + ^ + J^W + —y (52) 

where v' K is given by (1451) . Relations ( 1501) and ( 15T1) imply the condition determining 
extrema X r (b) of the radial epicyclic frequency profile: 

a ; = -^Q r (i,6). (53) 

In figure [U we show curves B^(x = X T ), k £ {1,2}, implicitly determined by 
relation ([53]) : index k denotes different branches of the solution of (!53|) : 

B k r = ^ (15 T V32x - 31) . (54) 
lb 

For all possible values of the brany parameter b < 1, the radial epicyclic frequency 
has one local maximum for braneworld RN neutron stars (black holes) that is always 
located above the marginally stable orbit x ms (see figured]). 

In the case of RN neutron stars with b > 1 (naked singularity type spacetimes), the 
situation is complicated, and the discussion can be separated into four parts according 



to the parameter b (see figures 1(b) and [3]): 
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(a) 1 < b < 5/4. 

The radial epicyclic frequency v x has one local maximum as usual in RN spacetimes. 

(b) 5/4 < b < b m = 1.313. 

The radial epicyclic frequency z/ r has one local maximum and one local minimum. 

(c) b m <b<b c = 1.381. 

The radial epicyclic frequency v x has no local extrema (for all b > b m ) since the 
condition x > b has to be satisfied, and for these values of b there is X r < x = b. 

(d) b > b c . 

The maximum of the Keplerian frequency profile could be relevant, since it can 
appear outside the internal uniform density configuration given by be (x). 

We can summarize that in braneworld RN spacetimes: 

• the Keplerian frequency is a monotonically decreasing function of radius for the 
whole range of the brany tidal charge parameter b in astrophysically relevant radii 
above the photon circular orbit (or the neutron star surface); an exception exists 
at spacetimes with b > 1.381, where the photon circular orbit is not defined. 

• the radial epicyclic frequency has a local maximum for RN black hole type 
spacetimes with b < 1, and vanishes at the innermost stable circular geodesic; 

• for RN naked singularity type spacetimes the behaviour of the radial frequency is 
different; a detailed analysis shows that the radial epicyclic frequency can have one 
local maximum for 1 < b < 5/4, two local extrema for 5/4 < b < 1.313 or no local 
extrema for b > 1.313. 

4-2. Orbital models of QPOs 

A number of low-mass X-ray binaries containing a neutron star show quasiperiodic 
oscillations (QPOs) in their X-ray flux, i.e., peaks in the Fourier variability power 
density spectra (PDS). Frequencies of some QPOs are in the kHz range corresponding to 
frequencies of the orbital motion close to the central neutron star. Such so-called high- 
frequency (kHz) QPOs span large frequency range of 200 - 1400 Hz. The two distinct 
modes following their own correlation between frequency and properties (amplitude and 
quality factor) of the peak are observed in this range (see, e.g., [2H EH [25] [26], [55], 156]). 
They are called lower and upper QPO because the frequencies of upper QPO, z/u, 
is higher than the frequency of the lower QPO, z/ L , when both modes are detected 
simultaneously; we call them twin peak QPOs in such situations [57] [36] [33] [58] . It was 
shown that the twin peak QPOs are clustered around frequencies corresponding to the 
frequency ratio of small integers [32] [33]. In some atoll sources, namely 4U 1636—53 
and 4U 1728, the frequency ratios are clustered around 3:2 and 5:4, in others (4U 1608 
and 4U 0614) they cluster around 3 : 2 only, and in 4U 1820 and 4U 1735 they cluster 
around ratio 4:3 [37]. In the best known Z-source Sco X-l the frequency ratios cluster 
around 3 : 2 and 5 : 4 [59] and in the exceptional Z-source Circinus X-l, the frequency 
ratio is clustered around 3:1 at relatively low frequencies ~ 100 Hz [60J. 




Figure 3. The Keplerian and radial epicyclic frequency profiles for some representative 
values of the brany parameter b in braneworld RN black hole (neutron star) (a) -(c), 
the extreme braneworld RN black hole case (d), RN naked singularity spacetimes (e) - 
(1). The radius of the outer event horizon Xh, the limiting photon orbit x p h and the 
marginally stable x ms circular orbits are denoted by a dotted vertical line (if these 
orbits exist at the given spacetime). Note that the condition x > b has to be satisfied, 
the dashed line in (e) - (1) represents the radius where x = b. The grey solid line 
denotes the critical radius Xq representing the limit on the radius of neutron stars. 
So the grey area is astrophysically irrelevant. (1) is the special case with z^k = v x at 
x = Xk = 2 for b = b SI = 1.5. 
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Several models have been outlined to explain observational data of the neutron 
star kHz QPOs assuming that their origin is related to orbital motion near the inner 
edge of accretion discs around the neutron stars. The RP model of Stella and Vietri 
[35] introduces the QPOs representing a direct manifestation of modes of a relativistic 
epicyclic motion of radiating blobs in the innermost parts of the accretion disc. In this 
model, the upper and lower frequencies are identified in the following way: 

u v = z/ K , v h = z/ K - z/ r . (55) 

The resonance model introduced by Abramowicz and Kluzniak [HH E2] (and earlier 
in other connections by Aliev and Galtsov [50]) assumes a non-linear resonance (forced 
or parametric) of accretion disc oscillations with vertical and radial epicyclic frequencies 
vq and u r . However, in spherically symmetric spacetimes assumed here u$ — z/k, and we 
identify 

z/rj = z/ K , vl = v x . (56) 

The observed frequency ratio z/u/^l decreases with increasing frequencies as x — > 
x ms , but u K /u r increases as x — > x ms (see [58] for details). On the other hand, 
frequency relations implied by the RP model yield a trend, which is in good accord 
with observation. In other words, frequency relation of the RP model satisfies in a 
natural way the general tendency of twin peak kHz QPOs as the ratio of z/k : {vk — Vt) 
decreases with increasing frequencies z/k and z/k — ^ r for x approaching the inner edge 
of the disc at x — x ms , where z/ r — > and z/k : (z/k — z/ r ) — ► 1. The clustering of the 
observational data could then be explained by non-linear resonant phenomena between 
the Keplerian and periastron oscillatory motion in the framework of the multiresonant 
model [58j H3]. Note, however, that this frequency relation can coincide with those 
implied under consideration of a forced resonance between the epicyclic frequencies 
with observed higher eigenfrequency z/u = z/k = ug and lower combinational frequency 
— v k — v T . In the spherically symmetric spacetimes, they coincide with radial m = 1 
and vertical m = 2 disc oscillation modes. Here, we focus for simplicity on the RP 
model. 



4-3. Resonance radii 

Radial profiles of the Keplerian and radial epicyclic frequency in the external spacetimes 
of the black hole type (b < 1) enable us to find the observed resonant points (3:2, 4:3, 
5 : 4) and relations for the resonant radii in the framework of both the relativistic 
precession model frequency relations z/k : (^k — v*) arid the epicyclic frequency relations 
u^:u r . The same relations hold in the case of external spacetimes of the naked singularity 
type with 1 < b < 5/4, where the behaviour of the Keplerian and radial epicyclic 
frequencies is the same as in the black hole type spacetimesjjj] 

§ We do not consider the internal branch of stable orbits since it is too close to the limit given by be- 
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Figure 4. The resonant radii for ratios of small integers for the RP model : (Vk — v-i) 
(a) and for the model of direct resonances j^k : v T (b). The grey lines represent the 
frequency ratios v\j : = 3:1, 2:1, 5:3, 3:2 (dotted line), 4:3 and 5:4 subsequently 
right to left in (a), left to right in (b). For 6 = and b = —5/4 the resonance radius 
where iaj ■ = 2 : 1 coincides with the radius of the radial epicyclic frequency local 
maximum (at x = 8 for b = 0, at x = 10 for b = —5/4). 
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The resonant radii x n:m , where the frequency ratio : (z/k — v r) given by the RP 
model takes the value n:m = p, are implicitly given by the condition 

b = b K ^\x, P ) EE 6 RP (x,p) 

(p - l) 2 x 2 + 9p 2 x =F xJx 2 (p - l) 4 + 2p 2 x (9p 2 - 2p + 1) - 15p 4 . 

S V • (57) 

For the model of direct resonances of the epicyclic frequencies we arrive at the relation 

b = b K/r (x, p\ ee ——7: \9p 2 + xt ^x 2 + 2p 2 x + p 4 (16x - 15)] . (58) 
8p z L 

These models are valid at x > x ms , but their behaviour is inverse as x — > x ms , since 
z/ K : (^k — ^r) —* 1 there, while z/ K : v T diverges there since v v {x = x ms ) = 0. 

The resonant radii for ratios of small integers allowing for efficient resonant 
phenomena [16] are given in figure HI Note that b^/ 2 (x) = b^.{ r (x), b^.{ ~ T \x) = 
bf/i '{ x )> ^h!i T \ x ) — bf/ 2 r {x). We stress that to model the data of twin peak QPOs 
in neutron star systems, we can use the RP model, but the direct epicyclic frequency 
model gives trends just opposite to what is observed. 

Using the RP model, we can conclude that the resonant points are located in the 
region between the marginally stable orbit and the local extremum of the epicyclic 
frequency, i.e., in the innermost part of the accretion disc. The frequency ratios ~ 3:1, 
observed in the Z-source Circinus X-l [60] could still be explained in the framework 
of the RP model, but the resonant points must be located above the maximum of the 
radial epicyclic frequency. The fitting of the observational data will be discussed in the 
following section, but for our rough estimates we shall concentrate on the behaviour 
of the frequency-frequency relations, postponing the detailed studies of the resonant 
phenomena in individual sources to future studies. 



4-3.1. Strong resonance radii. Here, we would like to mention some special resonance 
conditions that could be realized in the RN spacetimes of the naked singularity type. 

For RN naked singularity spacetime with b > 1.42, the Keplerian and radial 
epicyclic frequencies could be equal at a specific radius (see figure [3]). The condition 
VK.(x,b) = u T (x,b) implies the relation 

b = b ST (x) = |(9 + x =f y/x(x + 18) - 15) (59) 

determining in an implicit way the strong resonant radii, where the resonance between 
the oscillations with Keplerian and radial epicyclic frequencies could be strongest, 
and the frequency scatter allowed for the resonance is largest. The function b ST (x) 
is illustrated in figure [U Note, however, that the above relation is not relevant in case 
of neutron star systems. 



5. QPOs and RP model testing the braneworld models 

We have summarized the twin peak QPOs data for the atoll sources and Z-sources Sco 
X-l and Circinus X-l in figure [5] using the data accumulated in [25]. Clearly, there is 
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a general tendency in the vy-vi, frequency relations that can be immediately deduced 
from the observational data: i^j/^l decreases with (^u) increasing as stressed in the 
previous section. 

We shall use the most frequent RP model in order to obtain rough restrictions 
on the allowed values of the tidal charge (and brane tension) related to neutron star 
systems giving the observational data demonstrated in figure El 

Assuming that the trend observed in the frequency relations passing through the 
clusters of the neutron star twin peak kHz QPOs is related to the modes of the RP 
model, we make two kinds of fitting the data by the brany neutron star models. From 
these fittings we deduce rough constraints to the tidal charge and brane tension related 
to the neutron star systems generating the observed data. 

First, we assume negative values of the brany tidal charge (and positive brane 
tension A > 0) and consider any RP model with values of b < touching any 
observational point. Thus, we exclude the values of b < that correspond to the RP 
frequency relations located outside the observational data (see figure [5]). The results 



are characterized by M = 1.6 M Q and b £ (—2.8,0) (see figure 5(a)). Assuming the 
canonical neutron star mass M = 1.4 M Q , we obtain the acceptable interval of tidal 
charge b G (—3.6, 0) as demonstrated in figure 5(b)| 



Second, we allow both positive and negative values of b (and brane tension A). Such 
a case we represent by an estimate done in a symmetric way, i.e., we find such a neutron 
star mass that the observational data are scanned for this mass parameter by RP fit 
relations with brany tidal charge of the same size in both positive and negative values; 
we have found M = 2M Q and 6 G (—1.2, +1.2) as demonstrated in figure [6lfjj 1 

The presented results are used to obtain estimates on restrictions of the tidal charge 
B using the relation 

/ GM (i) \ 2 , s 

\B\<\B®\ = \H\\-^) ■ (60) 

We then obtain the following restrictions: 

(a) Assuming negative tidal charge with |6(i)| = 2.8 and M(i) = 1.6 M , we obtain 

\B\ < \B m \ = 1.562 x 10 n cm 2 , (61) 
while the canonical mass M( 2 ) = 1.4 M and |6(2)| = 3.6 yield 

\B\ < (5(2)1 = 1.537 x 10 11 cm 2 . (62) 

(b) Allowing both positive and negative values of the tidal charge, we obtain for 
M(3) = 2M Q and \brg)\ = 1.2 the estimate 

\B\ < 15(3)1 = 1.046 x 10 n cm 2 (63) 

|| Note that the range of b < 5/4 enables us to limit our consideration to the RN naked singularity 
type spacetimes with the same behaviour of the Keplerian and radial epicyclic frequencies as in the 
black hole spacetimes; for example, the strong resonance radii are irrelevant. 
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that is naturally lower in comparison with the previous two cases because of the 
scatter of b into both positive and negative values that reduces the magnitude of 
the dimensionless tidal charge by a factor of ~ 1.5. 

The brane tension can be estimated using the relation (J9]), which yields 

\M>\h)\-— (64) 

that can be transformed to the relation 

l A l > \ X (i)\ = \ h I I ~ n—\ X (i)- ( & 5) 
\0(i)\ r G \b(i)\ 

Using the limit on compactness of the uniform configuration R/tq > i?c/ r G> where 
Rc(b)/rc = Xc{b) is implicitly given as a function of b by relation fl36|) . we arrive at 
the general estimate on the brane tension: 

1^1 > 1^(0 1 = TI r " = T^ x c(%))- (66) 

Hi) I r G |0(»)| 

For the canonical values of the neutron stars, namely M = 1.4 M , R = 10 km and 
p ~ 6.6 x 10 14 g/cm 3 , we find 

|A| > |A can | = 2.7x 10 15 g/cm 3 . (67) 

In order to make a more precise restriction of A, we have to estimate the dimensionless 
surface radius X of the configuration for given bu\ and Muy For bu\ > —1.5, the estimate 
should be related to bc(x), for b^ < —1.5 we have to use b^{x). 
We find 

M \ 2 f 1 Y , jrjrj „ 1n l7. 



W>M -lM\j%) KWJ """"^cm'. (68) 

For Mm = 1.6 M , |6m| = 2.8 and choosing Xm = 4, we find 

| Ai| = 3.9 x 10 15 g/cm 3 , (69) 

which is comparable to our simple canonical estimate (1671) . For M($\ = 2 M , |6( 3 ) | = 1.2 
and choosing X( 3 ) = 3, we find 

|A 3 | = 1 x 10 16 g/cm 3 . (70) 

This estimate exceeds the canonical one by a factor of ~ 3.7 because of the lowering of 
the dimensionless tidal charge limit. 



6. Conclusions 

Using the RP model of kHz QPOs observed in neutron star X-ray binaries modified by 
the hypothetical tidal charge of the neutron stars as implied by the braneworld uniform 
star model, we are able to put rough restriction on the tidal charge magnitude and brane 
tension magnitude related to the considered group of atoll sources and some Z-sources. 
This is the first attempt at an estimate coming from the consideration of effects in the 
extremely strong gravitational fields. 
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b e (-2.8; 0) 




Figure 5. The relativistic precession model (with negative tidal charge) fits for atoll 
sources, Sco X-l and Cir X-l for various values of the brany parameter 6; the curves are 
spaced by 0.2 in b. The grey dashed lines represent the frequency ratios v\j : = 3 : 1, 
2:1, 5:3, 3:2, 4:3 and 5 : 4 subsequently left to right. 




Figure 6. The relativistic precession model (with positive and negative tidal charge) 
fits for the same ensemble of data as in figure Note that the M = 2 M Q , b = RP 
frequency relation is commonly used for the fitting of observed QPO data (see, e.g., 

1251). 



The presented method is rough but robust, since the effects of neutron star rotation 
are not important. According to (j4Tj) . the effect of the tidal charge is given by the term 
~ b/r 2 . On the other hand, the leading term giving rotational effects is ~ a/r 3 / 2 
and, moreover, the limits on the rotational parameters taken from observed rotational 
frequencies of neutron stars [36] put a strong limit of a < 0.3, so that its influence 
onto the character of the fitting frequency-frequency curves is small as shown by Torok 
et al. [36]. Of course, fitting the observational data of individual sources by precise 
methods will need the rotational parameter for fine tuning!^ 

Assuming negative values of B and positive values of A, we obtained rough estimates 
\B\ < 1.5 x 10 11 cm 2 and |A| > 3 x 10 15 g/cm 3 for both M ~ 1.6 M and the canonical 
case of neutron stars with M ~ 1.4 M Q . On the other hand, from the limits allowing 
both positive and negative values of b (A), starting from the standard mass estimate 
M ~ 2M related to QPOs [21], we obtained a lower estimate of \B\ < 10 n cm 2 , but 
higher values of |A| > 10 16 g/cm 3 . 

These values should be confronted with estimates given by Bohmer et al. |12j . 
concerning the effects of tidal charge in the limit of weak gravitational field related to 



% The precise estimates of individual sources probably require also inclusion of the non-geodesic effects 
caused by the neutron star magnetic field [63 . 
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the solar system. They considered the perihelion precession (Mercury), deflection of 
light near the edge of Sun, and the radar echo delay experiment. Clearly, the most 
convenient for comparison seems to be the first phenomenon, being of the same nature 
as the epicyclic motion of blobs in the relativistic precession model. The weak field 
restrictions are [32] 

\B\ < 5 x 10 8 cm 2 , A > 7 x 10 13 g/cm 3 . (71) 

Note that the rough strong field regime model gives estimates of the tidal charge higher 
than the weak field regime in more than two orders of magnitude. However, as stressed 
in [12], a substantial part of the perihelion shift could correspond to the effect of solar 
oblateness, modifying thus, in principle, very strongly the estimates of B and A. On the 
other hand, the restrictions implied in [12] by the light deflection and radar echo delay 
are \B\ < 10 19 cm 2 and |A| > 10 17 g/cm 3 , which are substantially, by orders, higher than 
those given by rough estimates of QPOs observational data in X-ray binaries. We can 
thus conclude that the restrictions on \B\ and |A| implied by rough fitting of the neutron 
star X-ray binaries kHz QPOs data, i.e., in the strong gravitational field regime, give 
better fits as compared to those implied by the weak-field limit, if these results are not 
masked by the effects of the solar oblateness. Clearly, further and more detailed studies 
of all these phenomena are necessary. Moreover, there is a variety of other vacuum 
solutions of spherically symmetric static field on the brane [6U [65j EE] that deserve 
attention and have to be tested separately since the Weyl tensor in the complete 5D 
theory is unknown yet. 

We plan to make further, more sophisticated comparisons of the plausible kHz 
QPOs models for neutron star X-ray binaries, including the RP model and the resonance 
disc oscillation ones, and compare these models with data obtained for individual 
sources. Precise fitting of these individual source data to the models including the 
influence of the tidal charge could, in principle, bring some new interesting results and 
limits on the tidal charge and brane tension although we cannot expect substantial shift 
of the rough estimates presented here, but rather better fits of observational data and 
more realistic estimates of the neutron star mass, more close to the expected canonical 
value. In fact, generally, we can conclude that the presence of a negative tidal charge 
enables lowering of the neutron star mass as can be seen from our rough estimates. 
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